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ABSTRACT 



OO 



Context. The ratio of pulsation to radial velocity (the projection factor) is currently limiting the accuracy of the Baade- 
Wesselink method, and in particular of its interferometric version recently applied to several nearby Cepheids. 
Aims. This work aims at establishing a link between the line asymmetry evolution over the Cepheids' pulsation cycles 
and their projection factor, with the final objective to improve the accuracy of the Baade-Wesselink method for distance 
q , determinations. 

Methods. We present HARPS* high spectral resolution observations (R — 120000) of nine galactic Cepheids : R Tra, 
S Cm, Y Sgr, j3 Dor, £ Gem, Y Oph, RZ Vel, £ Car and RS Pup, having a good period sampling (P — 3.39d to 
P = 41.52d). We fit spectral line profiles by an asymmetric bi-gaussian to derive radial velocity, Full- Width at Half- 
Maximum in the line (FWHM) and line asymmetry for all stars. We then extract correlations curves between radial 
velocity and asymmetry. A geometric model providing synthetic spectral lines, including limb-darkening, a constant 
FWHM (hereafter ac) and the rotation velocity is used to interpret these correlations curves. 

Results. For all stars, comparison between observations and modelling is satisfactory, and we were able to determine the 
projected rotation velocities and ac for all stars. We also find a correlation between the rotation velocity (V ro t sini) and 
the period of the star: V ro t sini = (—11.5 ±0.9) log(P) + (19.8 ± 1.0) [km. s -1 ]. Moreover, we observe a systematic shift in 
observational asymmetry curves (noted 70), related to the period of the star, which is not explained by our static model 
! 7o = ( — 10.7 ± 0.1) log(P) + (9.7 ± 0.2) [in %] . For long-period Cepheids, in which velocity gradients, compression or 
shock waves seem to be large compared to short- or medium-period Cepheids we observe indeed a greater systematic 
shift in asymmetry curves. 

Conclusions. This new way of studying line asymmetry seems to be very promising for a better understanding of Cepheids 
atmosphere and to determine, for each star, a dynamic projection factor. 
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1. Introduction The spectral line profile, in particular its asymmetry, is 

critically affected by the dynamical structure of Cepheids' 

Long-baseline interferometers currently provide a new atmosphere : photospheric pulsation velocity (hereafter 

quasi-geometric way to calibrate the Cepheid Period- v pu i s ), velocity gradients, limb-darkening, turbulence and 

Luminosity relation. Indeed, it is now possible to deter- rotation. Thus, radial velocities measured from line profiles, 

mine the distance of galactic Cepheids up to lkpc with the hereafter V rad , include the integration in two directions : 

Interferometric Baade-Wesselink method, her eafter IBW ver the surface, through limb-darkening, and over the ra- 

method (see fo r e.g . Sasselov & Karovska (| 1994ft and dius, through velocity gradients. All these phenomena, ex- 

Kervella et al. (|2004|), hereafter Paper I). Interferometric cep t the rotation, are supposed to vary with the pulsation 

measurements lead to angular diameter estimations over phase. However, they are currently merged in one specific 

the whole pulsation period, while the stellar radius varia- quantity, generally considered as constant with time: the 

tions can be deduced from the integration of the pulsation projection factor p, defined as V pu i s = pV ra d- 
velocity. The latter is linked to the observational velocity 

deduced from line profiles by the projection factor p. In this The interferometric definition of the projection factor 
method, angular and linear diameters have to correspond ig of crudal imp0 r ta nce in the IBW method, as it can in- 
to the same layer m the star to provide a correct estimate duce a biag of up tQ g% Qn the deriyed distance ( Nar detto 
of the distance. at aL (fjool), Merand et al. p55^) ). Otherwise, the limb- 

darkening is also required to derive a correct estimation of 

* High Accuracy Radial velocity Planetary Search project de- the angular diameter of the star. With the latest gener- 

veloped by the European Southern Observatory ation of long-baseline interferometers, studying its phase- 
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Table 1. Observed sample of Cepheids sorted by increasing period. 



Name 


HD 


P(a) 


To (a) 


Nb. 


Nb. of 


my (b) 






[days] 


[days] 


of spectra 


cycles 




R TrA 


135592 


3.38925 


2 451649.96 


14 


15 


6.66 


S Cm 


112044 


4.68976 


2 451 645.64 


12 


3 


6.60 


YSgr 


168608 


5.77338 


2 451650.92 


17 


10 


5.74 


(3 Dor 


37350 


9.84262 


2 451643.54 


49 


3 


3.75 


C Gem 


52973 


10.14960 


2 451641.78 


50 


3 


3.90 


Y Oph 


162714 


17.12520 


2 451653.32 


7 


4 


6.17 


RZ Vel 


73502 


20.40020 


2 451633.58 


10 


3 


7.08 


I Car 


84810 


35.551341 


2 452 290.4158 


118 


2 


3.74 


RS Pup 


68860 


41.51500 


2 451644.22 


15 


3 


7.03 



a For I Car, the reference Julian date (To) and the pulsatio n peri od (P) used to compute the phase are from Szabados i|1989f ). 

For others stars we used ephemeris from Berdnikov et al. (|200lh . 
b The visible magnitude (ray) is from Berdnikov et al. (|2000f ). 



depen dence is of crucial impo rtance (Marengo et al. (2002), 
(2003), Nardetto et al. (120061 )). 

Line asymmetry was fir st ob served for short-period 
cephe ids by Sasselov et al. (|1989t ). Then, Sasselov et al. 
(|1990T ) studied the impact of the asymmetry on radius and 
distances determinations. The link between line profiles 
asymmetry an d the projection factor has been studied by 
Albrow et al. (1994). Finally, an er ror an alysis of the IBW 
method is given in Marengo et al. (|2004[ ). 

We present here a new original study of the line asym- 
metry using the very high spectral resolution of HARPS 
(R = 120000). We have observed 9 galactic Cepheids with 
periods ranging from P—3.39d to P=41.52d Radial veloc- 
ity, full-width at half-maximum (hereafter FWHM) and line 
asymmetry are presented for all stars in Sect.[2j 

Section[3]deals with modelling and Sect.[2]with observa- 
tions interpretation. Through a geometric model different 
definitions of the projection factor are proposed and com- 
pared in order to find the best procedure. Then the model is 
used to interpret observational radial velocity and asymme- 
try correlation curves. A set of parameters is thus derived 
for all stars. Taking into account the whole sample of stars 
we discuss general properties and in particular the period- 
dependencies. 

2. HARPS observations 

2.1. Journal of observations 

HARPS is a spectrometer dedicated to the search for ex- 
trasolar planets by means of radial velocity measurements. 
It is installed at the Coude room of the 3.6 meter telescope 
at La Silla. The resolution is R = 120000 and the average 
Signal to Noise Ratio we obtain over all observations in 
the continuum (292 spectra) is 300 per pixel. The observed 
sample of Cepheids is presented in Table [TJ 

We have used the standard ESO/HARPS pipe- line re- 
duction package with a special attention for the normaliza- 
tion process. We have noted on metallic line profiles of all 
stars a good reproduction from cycle-to-cycle. Therefore, 
spectra for a given star have been recomposed into an 
unique cycle. 

Using Kurucz models (|1992f) we have identified about 
150 unblended spectral lines. This first study considers only 
the unblended metallic line Fel 6056.005 A. 



2.2. A new estimator of the radial velocity, FWHM and 
asymmetry: the bi-gaussian 

Several methods have been used to measure radial veloc- 
ities of Cepheids, each having advantages and drawbacks. 
Among these methods there is the line minimum (usually 
determined via a parabolic fit to a few pixels near the bot- 
tom of the line) a gaussian fit (obviously not adequate for 
asymmetric lines), the line centroid, determined from the 
integration of the line profile (requires high Signal /Noise 
ratio), and the line bisector where one measures the width 
of the line at one or several depths. Our bi-gaussian ap- 
proach combines advantages of methods useful for low S/N 
data while providing information usually associated with 
high resolution and high S/N data (asymmetry). 

Radial velocity, full width at half-maximum (FWHM) 
and asymmetry have been derived simultaneously apply- 
ing a classical x 2 minimization algorithm between the ob- 
served line profile (S(X)) and a modelled spectral line profile 
(/(A)). The corresponding reduced x 2 is : 

1 ^ (5(A t )-/(A t )) 2 

^ ^ (1) 

with N the number of pixel in the spectral line, v the 
number of degrees of freedom and a(\i) = SNR *f(Ai) is 
the statistical uncertainty associated to each pixel. SNR is 
the estimate of the Signal to Noise Ratio in the continuum. 

The analytic line profile is defined by : 

IW-i-D ^^-^ a^x. (2, 

and 

/(A) = 1 - Dexp ( 41n2(A-A m ) 2 \ 

M ' \ (FWHM (1 - A)) 2 ) w 

with four free parameters : 

— _D, the depth of the line. This quantity has no dimension. 

— A m , the wavelength associated to the minimum of the 
line (in A) . The corresponding radial velocity is noted 
RV m . 
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- FWHM is the Full- Width at Half-Maximum in the line, 
also in A. 

— A is the asymmetry as a percentage of the FWHM. 

The 4 In 2 factor is to obtain a correct definition of the 
FWHM. Forcing asymmetry to zero in this minimization 
process is equivalent to fitting a gaussian to the line profile. 
In this case we can derive another type of radial velocity 
noted RV g . 

There are different ways to define the line asymmetry 
(see e.g. Sasselov et al. 1990; Sabbey et al. 1995). The 
advantage of the bi-gaussian method is that it offers the 
possibility to derive statistical uncertainties directly from 
the minimization process. Moreover, all parameters (RV m , 
FWHM, D and A) are fitted simultaneously leading to a 
very consistent set of information. The largest reduced x 2 
we obtain with this method is of about 10 corresponding to 
a SNR of 438, but in most cases we have a reduced x 2 ~ 1 
or 2 corresponding to a SNR ranging from 75 to 350. That 
means that our analytic model is well suited to the data 
quality. We note also that the reduced X 2 is n °t sensitive 
to the spectral line resolution. 

As an example, Fig. [T] presents line profile variation for 
(3 Dor together with the analytic spectral line profile. We 
find that the asymmetry is insensitive to the choice of the 
continuum. However, this one has to be correctly defined 
to derive correct values of the FWHM and line depth D. 

Another radial velocity definition, the centroid velocity 
(RV C ) or, the first moment of the spectral line profile, has 
been estimated as: 



RVc. = 



L, c s{\)d\ 



(4) 



Tables [3l [4] and [5] present the resulting values of RV g , 
RV m , RV C , FWHM, D, A, SNR and X 2 red together with 
the corresponding uncertainties computed from the fitting 
method. 




6055,0 6055,5 6056,0 6056,5 

Heliocentric wavelength (A) 



6057,0 



6057,5 



Fig. 1. Spectral line evolution of (3 Dor together with the 
modelled bi-gaussian (bold). Line asymmetry is clear. The 
vertical line at the top corresponds to a differential flux of 
0.3. Pulsation phases are given on the right of each profile. 



2.3. Radial velocity 

As indicated in the previous section, we can derive three 
types of radial velocity : the velocity associated to the gaus- 
sian fit (RVg), the line minimum (RV m ) and the barycenter 
of the spectral line (RV C ). Figure[2]shows these radial veloc- 
ity curves obtained in the case of (3 Dor. Figure[3]represents 
for each star of our sample, the RV m variation (arbitrary 
shifted). The solid lines are the interpolated curves using a 
periodic cubic spline function. This function is calculated 
either directly on the observational points (e.g. (3 Dor) or 
using arbitrary pivot points (e.g. RZ Vel). In the latter 
classical minimization process between observations 
and the interpolated curve is used to optimize the position 
of the pivot points. All the interpolated curves presented 
in this study are derived using one of these two methods. 
The only exception is Y Oph (too few points) for which we 
performed a linear interpolation. 

2.4. The Full-Width at Half-Maximum in the line 

Figure 2] presents the FWHM curve as a function of phase 
for all stars. We note that the largest FWHM values are ob- 
tained for the maximum contraction velocities. RS Pup, the 
longest period Cepheid of our sample, seems to present an 
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0,0 0,2 0,4 0,6 0,8 1,0 1,2 

phase 

Fig. 2. (3 Dor radial velocities obtained with different 
method: RV m (points), RV g (squares), and RV C (crosses). 
Statistical uncertainties at ±lcr are indicated but too small 
to be visualized. We can therefore see the impact of the 
choice of the method in the case of a very asymmetric line 
(Fig. 
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0,0 0,5 1,0 1,5 2,0 

Phase 



Fig. 3. Radial velocity curves (RV m ). Curves have been 
arbitrarily shifted vertically. The horizontal lines are the 
zero velocity in the stellar rest frame. Largest velocities are 
for receding motion. 



important compression or shock wave signature. Figure [5] 
presents line profile variation for this star. Unfortunately 
the phase coverage is not very good, but we can clearly see 
a strong increase of the FWHM at <fi = 0.83. Such phe- 
nome non ha s been already detected in (3 Cepheids (Fokin 
ct al. (12004 )1. 

2.5. Asymmetry 

Figure [6] shows the asymmetry variation for all stars. 
Generally speaking, the shape of the aymmetry curve is 
similar to the shape of the velocity curve RV m . 

As already mentioned in Sect. 12.31 the radial velocity ac- 
cording to the choice of the method considered is sensitive 
to the line asymmetry. Figure [7] shows the correlation be- 
tween the differences of radial velocity (AV = RV m — RV g ) 
and the asymmetry of the line. We have only presented here 
the case of £ Car and RS Pup. Each star presents a similar 
behavior. A typical difference in velocity of about 4km. s -1 
can be obtained for an asymmetry of 40% in extreme cases 
(Y Sgr and R TrA). The relation between the radial veloc- 
ity difference and the asymmetry is certainly affected by 
star characteristics (rotation, FWHM, velocity gradients) 
present in the line asymmetry. In particular RS Pup signa- 
ture is certainly affected by strong velocity gradient effects. 
The fact that the RV m and RV g radial velocities present 
such differences as a function of the pulsation phase is an 
additional difficulty concerning an average projection factor 
and its time-dependence determination. With the centroid 
estimator of the radial velocity (RV g — RV C or RV m — RV C ) 
results are quite similar. 

In next sections, we summarize all observational results 
in correlation diagrams between radial velocity and asym- 
metry. These correlations are interpreted using the geomet- 
ric model in order to determine some physical parameters 
of our stars and to obtain information about dynamical 
effects in Cepheids atmosphere. 

3. A toy model 

We consider a limb-darkened pulsating star in rotation with 
an one-layer atmosphere. Our model has four parameters : 

— the limb-darkening of the star : we consider a linear law 
for the continuum-intensity profile of the star defined 
by I{cos(Q)) = 1 — ?iy + uy cos(0), where wy is_the limb- 
darkening of the star in V band (Claret et al. ( 2000)). Its 
value is about 0.7 for Cepheids. 9 is the angle between 
the normal of the star and the line-of-sight. 

— the projected rotation velocity V^ tsini, where i is the 
angle between the line-of-sight and the rotation axis (in 
km.s -1 ). 

— the pulsation velocity (in km.s -1 ). 

— the width of the spectral line (in A), hereafter named 
<7c ■ It is the FWHM of the line with no pulsation nor 
rotation velocities. It is supposed to be constant with 
the pulsation phase. 

The velocity field is a combination of pulsation and ro- 
tation velocities. Through the Doppler effect, this field can 
be transposed into wavelengths, and weighted by the sur- 
face brightness (limb-darkening) to obtain the weighting of 
the spectral line. We have then to convolve it with the in- 
trinsic profile to obtain the synthetic spectral line profile. 
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R TrA 



SCru 



YSgr 



fiDor 




YOph 




RZVel 



- % 



I Car 



RSPup 




0,0 



0,5 



1,0 
Phase 



1,5 



2,0 



Fig. 4. FWHM versus phase for all stars. Curves have 
been arbitrarily shifted vertically. The horizontal lines cor- 
respond to a zero FWHM. Note the particular case of RS 
Pup, which may present the signature of an important com- 
pression or shock wave. RS Pup has the longest period of 
our sample. 





6055,0 6055,5 6056,0 6056,5 6057,0 6057,5 
Heliocentric wavelength (A) 



6058,0 



Fig. 5. Fel 6056.005 A spectral line evolution of RS Pup. 
The vertical line at the top corresponds to a differential 
flux of 0.2. We note the broadening of the line at <f> = 0.83 
which could be the signature of a strong velocity gradient 
(compression or shock wave) . 



The weighting or the synthetic spectral line profile are pre- 
sented in different cases in Fig. [5] 

We now consider a pulsation velocity curve defined by: 



tpuls(<&) = VmaxCOs(27r0i) 



(5) 



with a typical value for the maximal pulsation velocity 
of V max = 30 km.s . This relation which is a poor ap- 
proximation of the pulsation velocity curve is only used for 
the projection factors determination (see below). It has no 
incidence on the results (see Sect 14. 1[) . From the synthetic 
spectral line profiles, we perform a bi-gaussian fit to de- 
rive the four parameters described in Sect. I2.2I : D, A m , 
FWHM and A. Then we derive the RV m , RV S , RV C ve- 
locities, and the corresponding radial velocity-asymmetry 
correlation curves (hereafter RV-A plot). In Fig. [TDl the 
RV-A plots are represented for different values of the ac 
and rotation parameters. The limb-darkening (considered 
as constant with the pulsation phase) has a very small ef- 
fect in the weighting of the line profile and thus practically 
no impact on the RV-A plot. Applying a classical minimiza- 
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Fig. 6. Asymmetry against phase for all stars. Curves have 
been arbitrarily shifted vertically. The horizontal lines cor- 
respond to an asymmetry of zero. 
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Fig. 7. Difference between the radial velocity obtained 
with the line minimum and the gaussian fit methods as a 
function of the asymmetry in the case of t Car and RS Pup. 
Statistical uncertainties are provided for each point. Arrows 
indicate the direction and the origin (j) = of the curves. 
These relations are not linear and certainly affected by star 
characteristics (rotation, FWHM, velocity gradients...). 



tion process between the pulsation and radial velocities, we 
have also derived for each set of parameters the correspond- 
ing constant projection factors : p m = V ^y° , Pg = V ^y B and 

f c ~ RV C ■ 

Firstly, we note that the ac of the line and the rota- 
tion have different effects on the slope and/or shape of the 
correlation curves. 

Secondly, correlation curves are slightly different from 
one definition of radial velocity to another. But the inter- 
esting point is that the RV C velocity does not depend of ac 
and/or rotation. This behavior is clearly seen on diagrams 
[T0b andfTUH : the centroid projection factor p c is constant 
with the ac and the rotation while the gaussian and the 
minimum projection factors, p g and p m , are varying. For 
the Cepheids of our sample the centroid projection factor 
ranges from p c = 1.40 (it v = 0.64; R TrA) to p c = 1.38 
(tiv = 0.75; £ Car), through the following relation : 



Pc 



-0.18u v + 1.52 



(6) 



This relation is an linear approximation from the ge- 
ometrical model (see Fig. Note that the geometrical 
model does not contain the physics of the pulsations, and 
thus the relation may not hold when instead of uy a more 
realistic limb-darkening (taking into account hydrodynamic 
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<7c = 0.1A J*f / 1 
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15 



-15 



-30 

-100%, 



-50% 



0% 

Asymmetry (%) 



50% 
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1 00% 0,0 
(d) 



Gc = 0.25 A 
uV=0.7 




Vrot sini=20km/s i 

ttf J / Vrot 
i *\ Vrot sini 


>ini=0km/s 
= 10km/s 



1 00% 



0,2 0,4 
Oc(A) 




5 10 15 

Rotation (km/s) 



Fig. 10. Results of the geometric model of pulsating star, (a-b) The radial velocity-asymmetry correlation curves for 
different o~c, with no rotation and no limb-darkening (uniform disk). Points, squares and crosses correspond respectively 
to the RV m , RV g and RV C radial velocities. For clarity RV g and RV m are represented only for oc = 0.1 A. The solid 
lines are the interpolated curves using a cubic spline function. The corresponding projection factors are represented on 
diagram (b). (c-d) Same plots but for different values of the rotation. The uq and the limb-darkening are respectively of 
0.25 A and 0.7. These RV-A plot are used to interpret HARPS observations. 



effects) is used. In particular, hydrodynamic effects can re- 
sult in a much larger limb-darkening, especially at the wave- 
lengths corresponding to spectral line (see e.g. Marengo et 
al. 2003). 



This behavior is of great importance in the context of 
the IBW method. Indeed, the community has often used 
the p 9 = 1.36 value of the projection factor (Burki et al. 
(1982)) using the gaussian method instead of the centroid 
metho d. As seen here, and already pointed out by Burki 
et al. (|1982f ). this estimator is biased by the rotation ve- 
locity, even if Cepheids are supposed to be slow rotators, 
and also by the oq . We thus recommend the centroid based 
methods (spectral observable and p-factor) for the analysis 
of Cepheid radial velocities. For the present work, we have 
therefore chosen the RV C definition of the radial velocity. 
Even though this requires substantial S/N, its advantages 
outweigh the drawback of spending more telescope time to 
acquire the data. 



4. Interpretation 

4.1. Methodology 

Modeling results obtained in the previous section are now 
helpful to elaborate a strategy in a comparison of observa- 
tions and models. 

Firstly, the effective temperature T c g and the surface 
gravity log g have been used to derive the intensity profile of 
stars through linear limb-darkening coefficients uy of Claret 
et al. (|2000l ) (see Tabled]). 

Secondly, we determine the projection factor p c us- 
ing Eq. [S] The pulsation velocity is then derived through 
Vpuis = PcRV c , where RV C is the observational radial veloc- 
ity corrected from the heliocentric velocity given in Table 
2. The pulsation velocity V pu i s and the projection factor 
p c (see Table 2) obtained are not physically realistic, be- 
cause our model does not include dynamical effects and in 
particular velocity gradients in the atmosphere, neverthe- 
less this procedure imposes the surimposition of observa- 
tional and modelled radial velocity curves RV C . Moreover, 
as a very good agreement is observed for each phase (bet- 



8 N. Nardetto et al.: High resolution spectroscopy for Cepheids distance determination 

Table 2. Optimized parameters obtained for each sample Cepheid through the confrontation of HARPS observations 
with our geometric model 



stars 



R TrA S Cm Y Sgr /3 Dor C Gem Y Oph RZ Vel 



Car 



RS Pup 



Period 

mean T e // a [K] 
mean \og(g) a 



4 



[kr 



3.38925 
6354 
2.0 
0.6371 
-13.2 



4.68976 
5995 
1.9 

0.6541 
-7.1 



5.77338 
5350 
1.0 

0.7194 
-2.5 



9.84262 
5490 

1.8 
0.6999 

7.4 



10.14960 
5727 
1.5 
0.6721 
6.9 



17.12520 
5907 
1.5 
0.6514 
-6.6 



20.40020 
5537 
1.5 
0.6970 
24.1 



oc WW] ' 
Kot sini (e) [km.s x ] 

P m — RV m 



35.551341 
5091 
1.5 
0.7541 
3.6 



41.51500 
5143 
0.4 
0.7121 
22.1 



Pg 

Pc - 



RVc 



0.29 
15 
1.13 
1.28 
1.40 



0.27 
10 
1.23 
1.31 
1.40 



0.27 
16 
1.10 
1.26 
1.39 



0.23 

6 
1.23 
1.32 
1.39 



0.23 

6 
1.23 
1.32 
1.40 



0.20 

4 
1.23 
1.33 
1.40 



0.23 

3 
1.26 
1.34 
1.39 



70 (e) [%\ 
7C (g) [%} 
■yo-c (h) 



0.25 

7 
1.23 
1.31 
1.38 



0.30 
< 1 
1.31 
1.36 
1.39 



3.3 
3.1 
0.2 



0.7 

4.3 
-3.6 



2.0 
0.4 
1.6 



0.2 
2.9 
-2.7 



-2.4 
0.5 
-2.9 



-3.2 
1.4 
-4.6 



-6.9 
1.2 
-8.2 



-6.5 
0.6 
-7.1 



T e ff[K] and log(p), deduced from Gieren et al. (T998) for R TrA, S Scu Y O ph and RZ Vel. For Y Sgr, (3 Dor, £ Gem, I Car, 
and RS Pup these quantities have taken from Cayrel de Strobel et al. f (|l997h . l|200lf )\ 
u v from Claret et al. (^OOPf ) 

v-y from Galactic Cepheid database (online : http://www.astro.utoronto.ca/DDO/research/cepheids) 
Uncertainty on ac is of about 0.02 A. 
Uncertainty on Kot sin i is of about 1km. s _1 . 

7o [%] is the averaged value of the observational asymmetry curves. The associated statistical uncertainties are of the order of 
0.3%. 

7c [%] is the averaged value of the computed asymmetry curves. 
7o-c [%] is the average value of the O-C asymmetry curve. 
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Fig. 8. The weighting or the synthetic spectral line profile 
in different cases, considering (a) the pulsation velocity, (b) 
the limb-darkening, (c) the rotation and, (d) an intrinsic 
width for the line (cc). 



ter than 1%), it validates the use of a constant projection 
factor (p c ). We find also that the poor description of the 
pulsation velocity (Eq. [5]) used to derive p c has no incidence 



1,41 



1,40 



Pc 1 ,39 



1,38 - 



1,37 




0,60 



0,80 



Fig. 9. The projection factor corresponding to the centroid 
velocity (p c ) as a function of the limb-darkening parameter 
(uy). Dots are the results from the toy model and the solid 
line corresponds to the linear approximation (x 2 — 10~ 5 ). 



on the resulting modelled RV C curve. By this procedure, we 
can thus concentrate only on the asymmetry, m aking the 
interpretation easier. Note that Nardetto et al. ([2004]) al- 
ready gave an indication of the impact of velocity gradients 
on the projection factor, and thus on the distance deter- 
mination, in the case of 5 Cep (about —6%). In Table 2, 
we also indicate for each star the corresponding projection 
factors p g and p m for comparison. 
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Thirdly, o~c and V I0 t sin i are determined together from 
the observational RV-A and FWHM curves. We first con- 
sider the minimum of the observational FWHM curve to 
obtain an indication on the value of ac ■ We then find the 
rotation which gives the best slope and shape for the RV- 
A curve. But as the rotation has also an impact on the 
FWHM (about 0.02A), we have then to slightly readjust 
ac accordingly. By this process we finally find the best and 
unique values for ac and V ro t sini. 

The uncertainties on V ro t sin i and ac , associated to the 
minimization process, were estimated to be respectively 
lkm. s _1 and 0.02A. Similar uncertainties are found if one 
considers several metallic lines. Note however that our toy 
model is too simple to provide secure and precise values of 
the rotation, which is the most interesting parameter. In 
particular the broadening of the spectral line due to the 
macro-turbulence can ce rtainly affect our rotation values 
(Bersier & Burki (|1996D ). Nevertheless our principal and 
first objective is to probe the dynamical effects by a direct 
comparison of our static model with observations. 

4.2. Observations Versus Modelisation 
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Fig. 12. Same as Fig.[TT]but for RS Pup. RS Pup seems to 
be a non-rotating star as requested by the shape of its RV-A 
curve. Note also atypical points in observational RV-A plot, 
which can certainly be interpreted through the presence of a 
strong compression or shock wave in the stellar atmosphere. 



We now apply our methodology to each Cepheid of our 
sample. Results are indicated in Table [2l RV-A plot are 
represented on Fig. QT] and [T2] Note that RV-A plot de- 
duced from the model have been shifted in asymmetry to 
match the observations (this point is discussed in next sec- 
tion). For R TrA and Y Sgr, we can notice a very small 
slope for the RV-A plot and a very large value for the ob- 
servational FWHM. It indicates a large rotational velocity 
V TOt sini and a properly small value for ac (see Fig. ITUac). 
Thus, the corresponding gaussian and minimum projection 
factors (p g and p m ) are lower than for others stars (see Fig. 
[TUbd). Conversely, for Y Oph and RZ Vel the RV-A plot 
have relatively large slope while the observational FWHM 
is typical (about 0.3). This has a direct consequence on 
the rotation, which is then very small, and on the projec- 
tion factors (p g and p m ) which are then relatively large. 
Comparatively, S Cru, (3 Dor and £ Gem can be considered 
as intermediate cases. For £ Car and RS Pup, we obtain an 
atypical RV-A plot which is greatly shifted in asymmetry. 
For RS Pup, we obtain a specific RV-A plot characterized 
by a strong curvature which can be interpreted by our geo- 
metric model as a very slow rotation velocity V ro tsini < 1 
km.s -1 . Note that atypical points which are observed at the 
top of the RV-A plot are certainly due to dynamical effects 
since they corresponds to phases of outwards acceleration. 

4.3. Discussion 

As observed in the particular case of i Car and RS Pup, 
an important systematic shift in asymmetry can be present 
between observations and models. We define respectively 
7o and 7c the averaged value of the observational and 
computed asymmetry curves [in %}. Note that the phases 
are sampled in the same way for data and model. Results 
are indicated in Tabled We have also calculated for each 
star the residuals between the observational and computed 
asymmetry curves, noted O-C curves (Fig. [T31 . We define 
J0-C1 the average value of these residual curves. These O- 
C asymmetry curves contain the whole dynamical informa- 
tion present in the observational asymmetry, mainly : the 



limb-darkening variation in the spectral line and with the 
pulsation phase, the micro- and macro- turbulence, velocity 
gradient and temperature effects. For R TrA, S Cru, Y Sgr, 
RZ Vel and RS Pup, we note a bump in the O-C asym- 
metry curves which is approximately linked to the cross of 
the compression wave just after the maximum contraction 
velocity (see Fig. However (3 Dor, £ Gem and I Car do 
not present such bump, which may be interpreted as the 
presence of a very small compression wave. In the case of 
Y Oph the phase sampling seems insufficient to conclude. 
Consistent hydrodynamical model would be helpful to con- 
firm these results. 

7o, 7c and 70 -c are represented as a function of the 
pulsation period on Fig. 114b . The open squares represent 
7c- We want to emphasize here that our model produces 
asymmetry curves with non-zero average value. Indeed, it 
is a natural consequence of the shape of the observational 
radial velocity curve used to derive the pulsation velocity. 
We find a similar behavior for all stars independently of the 
period. 

The shifts obtained on the observational asymmetry 
curves (70) show a very interesting linear dependence with 
the logarithm of the pulsation period: 



7o = (-10.7 ± 0.1) log(P) + (9.7 ± 0.2) [in %] (7) 

Moreover we note that the dependence of 70-c with 
the pulsation period is very similar to the one of 70- We 
can conclude that this behavior is related to the dynamical 
effects in the atmosphere, which are not taken into account 
in our toy model. This can be explained by the fact that 
long-period Cepheids have extended atmosphere and con- 
sequently strong velocity gradient (see for example the case 
of RS Pup mentioned above). Thus, the line forming region 
can be seriously perturbed leadin g to a systematic shift 
in asymmetry (Albrow & Cottrell (1994)). However, such 
an interpretation remains tricky and needs confirmation. 
Forthcoming hydrodynamical models are likely to bring out 
important insight in this field. 
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Fig. 11. Radial velocity (RV C ) - asymmetry correlation curves for R TrA, S Cru, Y Sgr, [3 Dor, Q Gem, Y Oph, RZ Vel 
and I Car. Dots and bold curves correspond respectively to observations and models. The statistical uncertainties are 
indicated. Note that RV-A plot deduced from the model have been shifted in asymmetry. The small plot on each diagram 
correspond to the comparison of the observational (dots) and model (bold curve) FWHM. 
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Fig. 13. Difference of the Observational and Computed 
asymmetry curves (O-C curves) for each stars. Curves are 
arbitrarily shifted. The horizontal dotted lines corresponds 
to a zero asymmetry for each star. 



From results of Tabled it appears also that the pro- 
jected rotational velocity varies as a function of the pulsa- 
tion period (Fig.lHb). We obtain the following relationship: 



V Iot sini = (-11.5 ± 0.9) log(P) + (19.8 ± 1.0) [in km.s" 1 ] 

(8) 

The projected rotation is an important parameter which 
can be used, for example, to study evolution of Cepheids 
together with their mass loss. However, note again that our 
toy model does not include the physics of the pulsations 
and it is also very difficult to separate the rotation and 
macroturbulence effects in the resulting broadening of the 
spectral line. Thus this relation has to be considered very 
carefully as it is certainly model dependent. 



5. Conclusion 

We have presented HARPS high spectral resolution 
(R= 120000) observations of nine galactic Cepheids hav- 
ing a good period sampling (P = 3.39d to P = 41.52d). 
We fit spectral line profile with an asymmetric bi-gaussian 
to derive radial velocity, FWHM and line asymmetry for 
all stars. The presence of a very important compression 
or shock wave in the case of RS Pup, the longest pe- 
riod Cepheid of our sample has been identified. We have 
also translated the measured spectroscopic quantities into 
meaningful correlation curves between radial velocity and 
asymmetry. 

A simple geometric model providing synthetic spectral 
lines, including limb-darkening, the ere and the projected 
rotation velocity is then used to interpret these correlations 
curves. 

Firstly, we find that the centroid projection factor (p c ) is 
independent of o~c and the rotation velocity. This projection 
factor is thus certainly the best one to use in the context 
of the Baade-Wesselink method. 

Secondly, we find for each stars an optimized set of 
parameters which allows to reproduce observational ra- 
dial velocity - asymmetry correlation curves. In particu- 
lar, we find a dependence of the derived projected rota- 
tion velocities with the period of the star : V ro t sin i = 
(-11.5 ± 0.9) log(P) + (19.8 ± 1.0) [in km.s" 1 ]. 
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Finally, by comparing the outputs of our static models '/? Dor' on page[l| 
and the observed quantities, we gain access to dynamical '£ Gem' on pageQ] 
effects. In particular, we found that long-period Cepheids 'Y Oph' on pageQ] 
with strong velocity gradient, like RS Pup, have a system- 'RZ VeF on pageQ] 
atic shift in their asymmetry curve. We thus derived a linear l £ Car' on page[j] 
relation between the observational shift in asymmetry and 'RS Pup' on pageQ] 
the logarithm of the period : 70 = (— 10.7 ± O.l)log(P) + 
(9.7 ± 0.2) [in %] . A detailed interpretation of these em- 
pirical relation is very difficult, but forthcoming hydrody- 
namical models are likely to bring out important insight in 
this field. 

In conclusion, line asymmetry, which contains most of 
the physics involved in Cepheid atmosphere, is an impor- 
tant tool. But additional hydrodynamical considerations to- 
gether with a multi-lines study are now required to have a 
better understanding of the dynamical processes present in 
Cepheid atmosphere and in particular to determine realistic 
projection factors including velocity gradients. 

Acknowledgements. Based on observations collected at La Silla obser- 
vatory, Chile, in the framework of European Southern Observatory's 
programs 072.D-0419 and 073.D-0136. This research has made use of 
the SIMBAD and VIZIER databases at CDS, Strasbourg (France). 
We thanks David Chapeau for his helpful collaboration concerning 
computing aspects, Olivier Chesneau and Philippe Stee for their care- 
ful reading of the manuscript, as well as Vincent Coude du Foresto 
and Andrei Fokin for useful discussions. 



References 

Albrow, M. D., & Cottrell, P.L. 1994, MNRAS, 267, 548 
Berdnikov, L.N., Dambis, A.K., & Vozyakova, O.V. 2000, 
A&ASuppl. Ser. 143, 211 

Berdnikov, L.N., & Caldwell, J. A. R. 2001, ApJ, 7, 3 

Bersier, D., & Burki, G. 1996, A&A, 306, 417 

Burki, C, Mayor, M., & Benz, W. 1982, A&A, 109, 258 

Cayrel de Strobel, C, Soubiran, O, Friel E.D., et al. 1997, 

A&ASuppl. Ser., 124, 299 

Cayrel de Strobel, C, Soubiran, O, Ralite, N. A&A, 373, 159 (2001) 
Claret, A. 2000, A&A, 363, 1081 

Fokin, A., Mathias, Ph., Chapellier, E., et al. 2004, A&A, 426, 687 
Gieren, W. P., Fouque, P., & Gomez, M. 1998, ApJ,496,17 
Kervella P., Nardetto N., Bersier D., et al. 2004, A&A, 416, 941 
Kurucz, R. L. 1992, IAU Symp. 149: The Stellar Populations of 
Galaxies, 149, 225 

Marengo, M., Sasselov, D. D., Karovska, M., et al. 2002, ApJ, 567, 
1131 

Marengo, M., Karovska, M., Sasselov, D. D., et al. 2003, ApJ, 589, 
975 

Marengo, M., Karovska, M., Sasselov, D. D., et al. 2004, ApJ, 603, 
285 

Merand A., Kervella P., Coude du Foresto V., et al. A&A438 (2005) 
L9-L12 

Nardetto, N., Fokin, A., Mourard, D., et al. 2004, A&A, 428, 131 
Nardetto, N., Fokin, A., Mourard, D., et al. 2006, A&A, accepted 
for publication. 

Sabbey, C. N., Sasselov, D. D., Fieldus, M. S., et al. 1995, ApJ, 446, 
250 

Sasselov, D. D., Karovska, M. 1994, ApJ, 432, 367 

Sasselov, D. D., Lester, J. B., Fieldus, M. S. 1989, ApJ, 337,29 

Sasselov, D. D., Lester, J. B., Fieldus, M. S. 1990, ApJ, 362, 333 

Sasselov, D. D., & Karovska, M., 1994, ApJ, 432, 367 

Szabados, L. 1989, Communications of the Konkoly Observatory 

Hungary, 94, 1 



List of Objects 

'R Tra' on page I 
'S Cru' on page 1 
'Y Sgr' on pageQ 



N. Nardetto et al.: High resolution spectroscopy for Cepheids distance determination 
Table 3. HARPS observations results for R TrA, S Cm and Y Sgr. 



13 



JD C 


phase 


Cy. 


Sp. 






RV C 




FWHM 


D 


A 


SNR 


Xred 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 












( 


l) 


G) 


(k) 


CO 


















R TrA 






















202 


.53 





.09 


14 


1 


-26.74 ± 0.07 


-29.88 ± 0.17 


-25.43 ± 1.16 





.467 


± 





.009 





.14 


-30.3 ± 2.7 


231 


2.0 


206 


.53 





.27 


15 


1 


-17.70 ± 0.05 


-19.71 ± 0.13 


i a no _i_ n r\i 

-lb. 92 ± 0.91 





.452 


± 





.005 





.17 


-19.5 ± 1.5 


224 


1.6 


152 


.65 


0. 


.37 


1 


2 


-11.18 ± 0.03 


-10.41 ± 0.08 


-11.39 ± 0.56 





.449 


± 





.002 


0. 


.19 


7.4 ± 0.8 


241 


1.6 


203 


.55 


0. 


.39 


14 


1 


-10.08 ± 0.04 


-9.06 ± 0.11 


-10.47 ± 0.71 





.447 


± 





.003 


0. 


.18 


9.9 ± 1.1 


257 


1.1 


156 


.65 





.55 


2 


2 


-0.98 ± 0.04 


2.20 ± 0.09 


-2.15 ± 0.52 





.460 


± 





.005 


0. 


.19 


30.4 ± 1.3 


209 


1.0 


204 


.52 





.67 


14 


1 


3.29 ± 0.06 


6.70 ± 0.15 


2.02 ± 0.74 





.492 


± 





.008 


0. 


.17 


30.6 ± 2.0 


207 


1.4 


150 


.65 





.78 


1 


2 


2.26 ± 0.04 


5.89 ± 0.09 


0.89 ± 0.50 





.534 


± 





.005 


0. 


.16 


29.9 ± 1.1 


255 


1.3 


201 


.54 





.79 


14 


1 


1.34 ± 0.06 


4.64 ± 0.14 


0.30 ± 0.74 





.524 


± 





.007 


0. 


.16 


27.5 ± 1.6 


247 


1.7 


154 


.65 


0. 


.96 


2 


2 


-24.46 ± 0.05 


-26.46 ± 0.13 


-23.76 ± 0.79 





.481 


± 





.004 


0. 


.12 


-18.3 ± 1.4 


241 


1.2 


205 


.54 





.98 


15 


1 


-25.60 ± 0.08 


-28.18 ± 0.19 


-24.71 ± 1.24 





.503 


± 





.007 


0. 


.12 


-22.8 ± 2.1 


243 


1.4 


















S CV11 

O Vyl LI 






















zU/ 


.40 


u 




3 


1 


-21.58 ± 0.05 


-23.61 ± 0.13 


-20.47 ± 1.01 





.386 


± 





.005 


0. 


.16 


-23.6 ± 2.0 


230 


1.5 


151 


.56 





.11 


1 


1 


-19.32 ± 0.05 


-21.43 ± 0.12 


-18.20 ± 0.93 





■ O t o 


-L. 


u 


fin 6 ; 


n 
u 


. lo 


-25.6 ± 2.0 


214 


1.4 


203 


.49 


0. 


.18 


3 


2 


-15.82 ± 0.03 


-17.26 ± 0.07 


-14.98 ± 0.50 





.344 


± 





.002 


0. 


.21 


-19.0 ± 1.0 


224 


1.7 


156 


.63 


0. 


.19 


1 


1 


-15.65 ± 0.04 


-16.90 ± 0.10 


-15.07 ± 0.52 





.281 


± 





.004 


0. 


.19 


-19.8 ± 1.9 


221 


2.4 


152 


.63 





.34 


1 


1 


-6.94 ± 0.02 


-7.08 ± 0.06 


-6.54 ± 0.46 





.306 


± 





.001 


0. 


.26 


-2.1 ± 0.8 


255 


1.9 


153 


.57 





.54 


1 


1 


4.61 ± 0.03 


5.87 ± 0.08 


4.21 ± 0.62 





.351 


± 





.002 


0. 


.26 


16.5 ± 1.1 


209 


2.5 


205 


.47 


0. 


.60 


3 


1 


7.09 ± 0.03 


8.77 ± 0.06 


6.44 ± 0.59 





.386 


± 





.002 


0. 


.26 


19.9 ± 0.9 


269 


2.5 


154 


.64 


0. 


.76 


1 


1 


13.40 ± 0.05 


15.66 ± 0.12 


12.48 ± 1.06 





.454 


± 





.005 


0. 


.23 


22.5 ± 1.4 


181 


1.6 


206 


.48 


0. 


.82 


3 


1 


12.07 ± 0.04 


14.42 ± 0.09 


11.17 ± 0.59 





.469 


± 





.004 


0. 


.20 


22.4 ± 1.0 


285 


1.5 


150 


.63 


0. 


.91 


1 


1 


-9.32 ± 0.14 


-10.19 ± 0.35 


-8.58 ± 1.37 





.412 


± 





.008 


0. 


.17 


-10.1 ± 3.7 


87 


1.4 


202 


.49 


0. 


.97 


3 


1 


-19.91 ± 0.05 


-21.62 ± 0.12 


-19.05 ± 0.70 

Y Scrr 





.401 


± 





.004 


0. 


.15 


-19.2 ± 1.5 


287 


2.0 


204 


.63 


0. 


.12 


10 


2 


-16 53 ± 06 


-20 47 ± 15 


I Dgl 

-15 07 ± 86 





.485 


± 





.011 


0. 


.15 


-36 4 ± 3 

'J V / . I — 1 — 'J ■ V / 


160 


1.1 


152 


.80 


0. 


.14 


1 


2 


-15.07 ± 0.04 


-18.08 ± 0.09 


-14.02 ± 0.56 





.488 


± 





.004 


0. 


.16 


-27.2 ± 1.3 


251 


1.4 


205 


.67 





.30 


10 


1 


-6.93 ± 0.05 


-8.45 ± 0.12 


-6.56 ± 0.65 





.472 


± 





.003 


0. 


.19 


-13.7 ± 1.1 


244 


2.6 


149 


.80 





.62 


1 


2 


9.37 ± 0.04 


12.58 ± 0.09 


8.13 ± 0.62 





.467 


± 





.005 


0. 


.21 


30.5 ± 1.3 


178 


1.0 


202 


.65 


0. 


.77 


10 


2 


18.50 ± 0.04 


23.46 ± 0.09 


16.08 ± 0.71 





.565 


± 





.007 


0. 


.19 


39.5 ± 1.4 


231 


1.5 


150 


.79 


0. 


.79 


1 


2 


18.31 ± 0.03 


23.31 ± 0.08 


15.98 ± 0.61 





.581 


± 





.006 


0. 


.18 


38.5 ± 1.1 


270 


1.8 


156 


.83 


0. 


.84 


2 


2 


13.31 ± 0.04 


18.07 ± 0.11 


11.56 ± 0.63 





.626 


± 





.007 


0. 


.16 


33.2 ± 1.2 


255 


1.6 


203 


.65 


0. 


.95 


10 


2 


-14.27 ± 0.04 


-16.27 ± 0.11 


-13.47 ± 0.50 





.530 


± 





.003 


0. 


.13 


-16.5 ± 1.0 


288 


1.4 


151 


.75 


0. 


.96 


1 


2 


-15.22 ± 0.05 


-17.50 ± 0.12 


-14.51 ± 0.61 





.517 


± 





.004 


0. 


.13 


-19.1 ± 1.2 


254 


1.1 



(a) JD C , average Julian date of observation defined by JD C — JD — 2453000 [in days]. 

(b) phase, averaged pulsation phase of observation. For ephemeris see Table [T] 

(c) Cy., pulsating cycle of the star corresponding to observation. 

(d) Sp., number of spectra associated to observation. Results corresponding to these spectra are averaged. 

(e) RV g , gaussian fit radial velocity and the associated error barre [in km.s^ 1 ]. 

(f) RVm, minimum radial velocity derived from the bi-gaussian fit [in fcm.s -1 ]. 

(g) RV C , radial velocity corresponding to the first moment of the spectral line [in km.s -1 ]. 

(h) FWHM, Full- Width at Half-Maximum derived from the bi-gaussian fit [in Angstroms]. 

(i) D, line depth derived from the bi-gaussian fit [no dimension]. Errors bars are not indicated but of the order of 10 
(j) A, asymmetry derived from the bi-gaussian fit [in percentage]. 

(k) SNR, observational spectral line signal to noise ratio. 

(1) X?edi reduced % 2 factor corresponding to the bi-gaussian fit. 
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Table 4. HARPS observations results for (3 Dor, £ Gem, Y Oph, and RZ Vel. See Tabled] for legend. 



JD C 


phase 


Cy. 


S P . 


RV e 


RVm 


RV C 


FWHM 


D 




SNR 


Xred 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


00 


(i) 


(j) 


(k) 


(i) 



/3 Dor 



21.68 


0.02 


1 


4 


1.70 ± 0.01 


0.99 ± 0.02 


2.17 ± 0.14 


0. 


.286 


± 





.001 


0.23 


-11.3 zb 0.4 


345 


3.7 


31.64 


0.03 


2 


3 


1.35 ± 0.01 


0.68 ± 0.02 


1.64 ± 0.12 


0. 


.275 


± 





.001 


0.23 


-11.0 zb 0.4 


404 


2.4 


32.68 


0.14 


2 


3 


-5.16 ± 0.01 


-6.59 ± 0.03 


-4.33 ± 0.19 


0. 


.318 


± 


0. 


.001 


0.24 


-20.3 zb 0.6 


298 


2.1 


23.64 


0.22 


1 


4 


-0.73 ± 0.01 


-1.42 ± 0.01 


-0.19 ± 0.11 


0. 


.280 


± 


0. 


.001 


0.30 


-11.4 zb 0.2 


423 


4.9 


33.61 


0.23 


2 


3 


0.15 ± 0.01 


-0.49 ± 0.02 


0.71 ± 0.12 


0. 


.275 


± 


0. 


.001 


0.31 


-11.0 zb 0.3 


443 


7.9 




u.oo 








y.ui zl u.ui 


O qr> _i_ n no 
y.yu zn u.uz 


q RC 4- n 99 
y.uo ziz u.zz 


0. 


.253 


4- 
zn 


n 
u. 


nni 

UU1 


u.ou 


z ±. i zn u.o 


oou 


9 r i 


15.62 


0.40 


1 


3 


16.16 ± 0.01 


16.85 ± 0.02 


15.85 ± 0.28 


0. 


.261 


± 


0. 


.001 


0.35 


12.3 =b 0.4 


262 


2.4 


25.68 


0.42 


2 


3 


18.10 ± 0.01 


18.80 ± 0.02 


17.86 ± 0.24 


0. 


273 


± 


0. 


.001 


0.33 


11.8 ± 0.3 


399 


3.3 


35.64 


0.44 


3 


2 


19.09 ± 0.01 


19.87 ± 0.02 


18.80 ± 0.35 


0. 


.290 


± 


0. 


.001 


0.33 


12.3 zb 0.4 


337 


2.0 


16.67 


0.51 


1 


3 


24.95 ± 0.01 


26.53 ± 0.02 


24.41 ± 0.38 


0. 


.347 


± 


0. 


.001 


0.28 


20.5 ± 0.4 


352 


3.2 


26.59 


0.52 


2 


2 


25.48 ± 0.01 


27.09 ± 0.02 


24.91 ± 0.34 


0. 


.359 


± 


0. 


.001 


0.27 


20.2 zb 0.3 


473 


4.9 


36.64 


0.54 


3 


2 


26.61 ± 0.01 


28.52 ± 0.04 


25.83 ± 0.54 


0. 


.388 


± 


0. 


.001 


0.26 


22.1 ± 0.5 


336 


2.2 


17.69 


0.61 


1 


3 


27.57 ± 0.02 


30.11 ± 0.04 


26.60 ± 0.51 


0.457 


± 


0. 


.002 


0.22 


24.8 zb 0.5 


303 


2.2 


37.64 


0.64 


3 


2 


25.54 ± 0.02 


27.46 ± 0.04 


24.86 ± 0.45 


0.473 


± 


0. 


.001 


0.20 


18.0 ± 0.5 


409 


3.1 


28.67 


0.73 


2 


3 


11.09 ± 0.01 


11.09 ± 0.03 


11.15 ± 0.20 


0.401 


± 


0. 


.001 


0.21 


0.1 zb 0.3 


456 


2.7 


29.63 


0.83 


2 


4 


2.24 ± 0.01 


1.65 ± 0.02 


2.60 ± 0.12 


0. 


343 


± 


0. 


.001 


0.23 


-7.9 ± 0.2 


472 


6.5 


30.59 


0.92 


2 


3 


2.36 ± 0.01 


1.71 ± 0.02 


2.77 ± 0.12 


0. 


.286 


± 


0. 


.001 


0.24 


-10.5 zb 0.3 


455 


5.0 














C, Gem 




















32.70 


0.04 


2 


3 


-3.82 ± 0.02 


-5.04 ± 0.04 


-3.05 ± 0.33 


0. 


.313 


zb 


0. 


.001 


0.27 


-18.2 zb 0.7 


196 


2.9 


33.62 


0.14 


2 


3 


-4.67 ± 0.01 


-6.02 ± 0.02 


-3.55 ± 0.19 


0. 


.292 


zb 


0. 


.001 


0.31 


-21.9 zb 0.4 


330 


7.7 


23.65 


0.15 


1 


4 


-4.00 ± 0.01 


-4.89 ± 0.02 


-3.35 ± 0.13 


0. 


.276 


zb 


0. 


.001 


0.31 


-15.0 =b 0.3 


338 


4.2 


34.65 


0.23 


2 


3 


1.53 ± 0.01 


1.23 ± 0.02 


1.97 ± 0.15 


0. 


.248 


± 


0. 


.001 


0.37 


-5.7 =b 0.3 


334 


6.0 


35.65 


0.34 


2 


3 


9.44 ± 0.01 


9.73 ± 0.02 


9.41 ± 0.18 


0. 


.239 


± 


0. 


.001 


0.39 


5.8 zb 0.3 


299 


2.6 


25.69 


0.35 


I 


3 


1 86 + 01 


11 99 + 03 


1 84 + 35 


0. 


.256 




0. 


.001 


0.37 


6.6 zb 0.4 


195 


1.9 


15.71 


0.37 


1 


3 


12.29 ± 0.01 


12.81 ± 0.02 


12.12 ± 0.25 


0. 


.259 


± 


0. 


.001 


0.38 


9.4 zb 0.3 


253 


1.6 


36.66 


0.43 


2 


2 


17.10 ± 0.01 


17.87 ± 0.03 


16.92 ± 0.46 


0. 


.304 


± 


0. 


.001 


0.34 


11.7 ± 0.5 


255 


2.3 


26.60 


0.44 


2 


3 


17.69 ± 0.01 


18.53 ± 0.02 


17.48 ± 0.29 


0. 


.308 


± 


0. 


.001 


0.33 


12.4 zb 0.3 


353 


3.1 


16.69 


0.46 


1 


3 


19.06 ± 0.02 


20.11 ± 0.04 


18.71 ± 0.50 


0. 


.324 


± 


0. 


.001 


0.31 


14.6 ± 0.6 


189 


1.4 


37.66 


0.53 


3 


2 


21.81 ± 0.01 


23.27 ± 0.03 


21.40 ± 0.45 


0. 


.397 


± 


0. 


.001 


0.28 


16.4 zb 0.4 


341 


3.1 


17.70 


0.56 


1 


3 


21.79 ± 0.01 


23.25 ± 0.02 


21.35 ± 0.33 


0.439 


± 


0. 


.001 


0.27 


14.9 zb 0.2 


446 


3.1 


28.68 


0.62 


2 


2 


16.11 ± 0.02 


16.81 ± 0.06 


16.02 ± 0.59 


0.449 


zb 


0. 


.001 


0.24 


7.0 zb 0.6 


243 


1.4 


29.64 


0.74 


2 


2 


6.25 ± 0.03 


5.98 ± 0.07 


6.53 ± 0.53 


0. 


.372 


zb 


0. 


.001 


0.26 


-3.3 =b 0.8 


169 


1.7 


30.60 


0.84 


2 


3 


1.64 ± 0.01 


1.06 ± 0.02 


2.09 ± 0.16 


0. 


.321 


± 


0. 


.001 


0.28 


-8.5 =b 0.3 


407 


7.2 


31.64 


0.94 


2 


3 


0.49 ± 0.01 


-0.13 ± 0.02 


0.98 ± 0.15 


0. 


.282 


zb 


0.001 


0.28 


-10.2 zb 0.3 


372 


8.3 


21.70 


0.96 


1 


5 


0.18 ± 0.01 


-0.54 ± 0.03 


0.59 ± 0.15 


0. 


.271 


zb 


0.001 


0.27 


-12.1 zb 0.5 


229 


2.3 














Y Oph 




















216.75 


0.29 


4 


1 


-9.75 ± 0.02 


-9.83 ± 0.05 


-9.72 ± 0.46 


0. 


.205 


± 


0. 


.001 


0.32 


-1.9 zb 1.1 


189 


1.2 


201.63 


0.41 


4 


1 


-4.48 ± 0.03 


-4.31 ± 0.08 


-4.57 ± 0.43 


0. 


.202 


± 


0. 


.002 


0.31 


4.0 zb 1.8 


119 


1.4 


150.78 


0.44 


1 


1 


-3.15 ± 0.02 


-2.87 ± 0.04 


-3.37 ± 0.27 


0. 


.223 


± 


0. 


.001 


0.33 


5.9 zb 0.7 


262 


2.9 


203.65 


0.53 


4 


1 


0.57 ± 0.02 


1.06 ± 0.04 


0.20 ± 0.22 


0. 


.238 


± 





.001 


0.31 


9.7 zb 0.7 


296 


4.3 


152.80 


0.56 


1 


1 


1.37 ± 0.02 


1.80 ± 0.04 


1.09 ± 0.24 


0. 


.244 


± 





.001 


0.30 


8.3 zb 0.7 


297 


3.8 


154.75 


0.67 


1 


1 


2.89 ± 0.03 


3.30 ± 0.07 


2.47 ± 0.41 


0. 


.263 


zb 


0. 


.001 


0.27 


7.4 zb 1.1 


208 


3.6 


156.71 


0.79 


1 


1 


-1.85 ± 0.04 


-1.62 ± 0.10 


-2.13 ± 0.63 


0. 


.259 


zb 


0. 


.002 


0.26 


4.1 zb 1.6 


142 


2.5 














RZ Vel 




















204.44 


0.00 


3 


1 


13.05 ± 0.23 


12.50 ± 0.58 


14.39 ± 4.04 


0. 


.588 


zb 


0. 


.012 


0.14 


-4.4 zb 4.0 


76 


1.3 


205.44 


0.05 


3 


1 


-0.47 ± 0.07 


-3.76 ± 0.16 


0.69 ± 0.78 


0.457 


zb 


0. 


.009 


0.16 


-32.0 zb 2.7 


210 


1.5 


206.44 


0.10 


3 


1 


-1.90 ± 0.08 


-6.21 ± 0.17 


-0.30 ± 0.96 


0.459 


zb 





.016 


0.17 


-42.7 zb 4.9 


162 


2.5 


150.49 


0.36 


1 


1 


18.65 ± 0.01 


18.46 ± 0.03 


18.76 ± 0.38 


0.219 


zb 


0. 


.001 


0.38 


-3.9 =b 0.5 


309 


7.8 


152.51 


0.46 


1 


1 


28.82 ± 0.01 


28.83 ± 0.02 


28.61 ± 0.60 


0.231 


zb 


0. 


.001 


0.41 


0.8 =b 0.5 


205 


2.8 


154.50 


0.55 


1 


1 


45.64 ± 0.02 


47.66 ± 0.05 


44.86 ± 1.46 


0.369 


zb 


0. 


.002 


0.32 


24.7 zb 0.9 


239 


3.2 


156.49 


0.65 


1 


1 


43.14 ± 0.05 


43.91 ± 0.14 


42.87 ± 2.31 


0.526 


zb 


0. 


.003 


0.23 


6.4 zb 1.1 


178 


1.8 


201.44 


0.86 


3 


1 


39.12 ± 0.04 


39.87 ± 0.09 


38.93 ± 1.45 


0.411 


zb 


0. 


.002 


0.24 


8.1 =b 1.0 


224 


2.7 


202.45 


0.90 


3 


1 


29.54 ± 0.06 


29.84 ± 0.15 


29.26 ± 1.09 


0.483 


zb 


0. 


.003 


0.14 


2.8 zb 1.2 


271 


1.2 


203.44 


0.95 


3 


1 


13.05 ± 0.23 


12.50 ± 0.58 


14.39 ± 4.04 


0.588 


zb 


0. 


.012 


0.14 


-4.4 zb 4.0 


76 


1.3 
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Table 5. HARPS observations results for I Car and RS Pup. See Tabl^3]for legend. 
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JD C 


phase 


Cy. 


S P . 


RV S 


RV m 


RV C 


FWHM 


D 


A 


SNR 


Xred 


(a) 


(b) 


(c) 


(d) 


(e) 


(f) 


(g) 


(h) 


0) 


(j) 


00 


(i) 



I Car 



37.65 


0.02 


1 


7 


-13.40 ± 0.01 


-15.70 ± 0.02 


-12.19 ± 0.18 


0.416 


± 


0.001 


0.26 


-25.4 ± 0.3 


354 


3.5 


40.63 


0.10 


1 


5 


-12.24 ± 0.01 


-14.22 ± 0.02 


-11.10 ± 0.17 


0.362 


± 


0.001 


0.31 


-25.2 ± 0.3 


333 


3.6 


47.69 


0.30 


1 


2 


-1.57 ± 0.01 


-2.51 ± 0.02 


-0.61 ± 0.24 


0.274 


± 


0.001 


0.42 


-16.4 ± 0.4 


275 


7.5 


48.62 


0.33 


1 


2 


-0.01 ± 0.01 


-0.76 ± 0.01 


0.59 ± 0.10 


0.257 


± 


0.001 


0.40 


-13.7 ± 0.2 


438 


10.3 


49.67 


0.36 


1 


2 


1.79 ± 0.01 


1.09 ± 0.02 


2.34 ± 0.13 


0.260 


± 


0.001 


0.40 


-12.7 ± 0.3 


374 


6.0 


15.72 


0.40 


1 


3 


4.92 ± 0.00 


4.33 ± 0.01 


5.42 ± 0.12 


0.276 


± 


0.001 


0.40 


-10.2 ± 0.2 


405 


7.0 


51.68 


0.41 


2 


4 


5.23 ± 0.00 


4.62 ± 0.01 


5.75 ± 0.13 


0.274 


± 


0.001 


0. 


.39 


-10.4 ± 0.2 


352 


6.0 


16.69 


0.43 


1 


3 


6.64 ± 0.01 


6.13 ± 0.02 


7.15 ± 0.21 


0.289 


± 


0.001 


0. 


.39 


-8.3 ± 0.3 


293 


5.6 


52.64 


0.44 


2 


2 


6.93 ± 0.01 


6.30 ± 0.02 


7.42 ± 0.18 


0.285 


± 


0.001 


0. 


.39 


-10.3 ± 0.3 


376 


5.4 


17.71 


0.46 


1 


3 


8.46 ± 0.01 


8.11 ± 0.01 


8.86 ± 0.15 


0.300 


± 


0.001 


0. 


.38 


-5.5 ± 0.2 


444 


8.7 


53.69 


0.47 


2 


2 


8.70 ± 0.01 


8.31 ± 0.02 


9.00 ± 0.18 


0.288 


± 


0.001 


0. 


.37 


-6.3 ± 0.3 


390 


5.0 


54.67 


0.50 


2 


2 


10.36 ± 0.01 


10.18 ± 0.02 


10.52 ± 0.22 


0.296 


± 


0.001 


0. 


.35 


-2.9 ± 0.3 


340 


4.7 


55.70 


0.53 


2 


2 


11.99 ± 0.01 


12.04 ± 0.02 


12.08 ± 0.20 


0.316 


± 


0.001 


0. 


.35 


0.6 ± 0.2 


418 


4.5 


20.84 


0.54 


1 


3 


14.11 ± 0.01 


14.65 ± 0.03 


13.93 ± 0.35 


0.335 


± 


0.001 


o. 


.33 


7.4 ± 0.4 


218 


2.2 


56.70 


0.55 


2 


2 


13.46 ± 0.01 


13.75 ± 0.02 


13.42 ± 0.25 


0.331 


± 


0.001 


0. 


.34 


4.0 ± 0.3 


378 


4.2 


21.85 


0.57 


1 


3 


15.01 ± 0.01 


15.66 ± 0.02 


14.76 ± 0.22 


0.358 


± 


0.001 


0. 


.32 


8.3 ± 0.2 


379 


3.6 


57.70 


0.58 


2 


2 


14.82 ± 0.01 


15.35 ± 0.02 


14.67 ± 0.27 


0.347 


± 


0.001 


0. 


.34 


7.0 ± 0.3 


384 


3.7 


58.71 


0.61 


2 


2 


16.00 ± 0.01 


16.83 ± 0.02 


15.70 ± 0.27 


0.365 


± 


0.001 


0. 


.33 


10.5 ± 0.3 


412 


5.3 


23.66 


0.62 


1 


4 


16.98 ± 0.01 


17.92 ± 0.02 


16.55 ± 0.20 


0.388 


± 


0.001 


0. 


.31 


11.1 ± 0.2 


428 


3.5 


24.85 


0.66 




4 


1 7 94 + 01 


1 Q 1 5 + 09 


1 7 AO + D 99 


0.389 


-(- 


0.001 


0. 


.30 


14.0 ± 0.2 


374 


3.6 


25.87 


0.69 


1 


4 


18.52 ± 0.01 


19.63 ± 0.02 


18.00 ± 0.19 


0.404 


± 


0.001 


0. 


.30 


12.4 ± 0.2 


492 


4.2 


26.85 


0.71 


1 


5 


19.03 ± 0.01 


20.27 ± 0.02 


18.50 ± 0.24 


0.426 


± 


0.001 


0. 


.31 


13.4 ± 0.2 


436 


5.9 


28.69 


0.77 


1 


5 


20.12 ± 0.01 


21.26 ± 0.02 


19.50 ± 0.22 


0.426 


± 


0.001 


0. 


.29 


12.2 ± 0.2 


433 


3.1 


29.65 


0.79 


1 


5 


20.33 ± 0.01 


21.61 ± 0.02 


19.63 ± 0.23 


0.445 


± 


0.001 


0. 


.29 


13.2 ± 0.2 


439 


4.7 


30.80 


0.83 


1 


7 


19.76 ± 0.01 


20.83 ± 0.01 


19.05 ± 0.19 


0.450 


± 


0.001 





.28 


11.0 ± 0.1 


446 


6.6 


31.66 


0.85 


1 


5 


17.77 ± 0.01 


18.65 ± 0.02 


17.15 ± 0.24 


0.458 


± 


0.001 


0. 


.27 


8.9 ± 0.2 


391 


5.3 


32.72 


0.88 


1 


5 


11.81 ± 0.01 


12.17 ± 0.02 


11.51 ± 0.20 


0.481 


± 


0.001 


0. 


.24 


3.4 ± 0.2 


431 


3.8 


33.63 


0.91 


1 


5 


3.92 ± 0.01 


3.06 ± 0.03 


4.31 ± 0.20 


0.483 


± 


0.001 


0. 


.23 


-8.2 ± 0.2 


371 


4.1 


34.67 


0.93 


1 


5 


-4.93 ± 0.01 


-6.89 ± 0.02 


-3.99 ± 0.15 


0.460 


± 


0.001 


0. 


.24 


-19.5 ± 0.2 


486 


8.3 


35.66 


0.96 


1 


7 


-9.85 ± 0.01 


-11.93 ± 0.02 


-8.81 ± 0.14 


0.434 


± 


0.001 


0. 


.24 


-21.8 ± 0.2 


421 


6.0 


36.65 


0.99 


1 


6 


-12.39 ± 0.01 


-14.56 ± 0.02 


-11.29 ± 0.21 
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